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We investigated the current- voltage characteristic of InAs/GaAs quantum dot inter- 
mediate band solar cells (QD IBSCs) with different n-type doping density in the QD 
layer. The n-type doping evidently increases the open circuit voltage, meanwhile 
decreases the short circuit current density, and leads to the conversion efficiency ap- 
proaching that of the control solar cell, that is the major role of n-type doping is to 
suppress the effects of QDs on the current-voltage characteristic. Our model adopts 
practical parameters for simulation rather than those from detailed balanced method, 
so that the results in our simulation are not overestimated. 
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I. INTRODUCTION 



Intermediate band solar cells (IBSCs) are expected to be the third generation photovoltaic 
devices because of the potential for high energy conversion efficiency.- - - In addition to the 
absorption between the conduction band (CB) and the valence band (VB) the same as single 
junction solar cells, there are transitions from the VB to the IB and from the IB to the CB, 
it is expected that the IBSCs have the advantages of large increase of short circuit current 
density, small decrease of open circuit voltage, and enhanced conversion efficiency. 

Quantum dots (QDs) are proposed to form the IB,- - - and another implementation is 
employing alloys with intrinsic IB.->2 By now the experiments show that the efficiency of 
InAs/GaAs QD IBSCs are lower than control GaAs solar cells without QDs due to the 
reduction of open circuit voltage and small increase of short circuit current density.-^^ 

N-type doping in the i-region are raised to half-fill the IB so that the IB can provide 
strong absorption in transitions from the IB to both the VB and the CB.— Modulated 5- 
doping is generally carried out with sheet density equal to the QD density,—"— the doped 
layer is a certain distance away from the QD layer, and there is approximately one electron 
on the dot level per quantum dot because there are two states with different spin. Recently 
direct doping Si into QDs not only can half-fill the IB but also lead to enhancement of the 
photoluminescence intensity by a dedicated design doping in the QD's assembling stage,— 
then this method is applied to the IBSC and increase in photocurrent is observed.— 

In this paper, we employ the drift-diffusion model to investigate the effects of InAs/GaAs 
QDs on current- voltage characteristic in IBSCs by changing n-type doping density in the QD 
layer. The main difference between this modeling and previous work done on IBSCs, is that 
previously the current-voltage characteristics are mostly based on the detailed balanced 
method,-^ 1 ^ 1 ^ material parameters such as the absorption coefficients, carrier mobilities, 
etc. are not involved. Although some works adopt drift-diffusion model for the simulation 
of the device structure, the generation and recombination via the IB are still deduced from 
the detailed balanced method,— 1 ^ because of the large difference between the theoretical 
and experimental results, it is necessary to adopt practical parameters. 
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II. THE THEORY MODEL 



A. Drift-diffusion model 

In solar cells, the electron and hole current equations include both drift and diffusion 
components, 

dip dn 

J n = enfx n - eD n — , (1) 

T dip dp 

J P = -ep» P - - eD p -, (2) 

where e is the electron charge, n and p are the electron and hole concentration, fi n and \x v 
are the carrier mobilities, which decrease to small values as doping density increases, D n 
and D p are the carrier diffusion constants determined by the Einstein relation, D n = —u n 
and D p = —u p . 

Drift-diffusion model includes Poisson equation, electron and hole current continuity 
equations,— ~— when considering the role of QDs, the three equations can be written as 



d_ 

dx 



=e[p-n-fN I + N+-Nj}, 



(3) 



^ = e {G cv + Gj- R cv - Rsrh), (4) 
-j^- = e(Gcv + Gj — Rcv — Rsrh), (5) 

where e is the permittivity of the medium, ip is the electrostatic potential, Ni is the number 
of the IB states per unit volume, / is the electron occupation factor in the IB, which will 
be studied later, Np is the ionized donors and iV^" is the ionized acceptors. 

Gov — J &(X)Io(X)R(X) exp[— a(X)x]d\, is the generation rate of carriers from the VB 
to CB, where Io(X) is the intensity of the light at x = 0, a(A) is the absorption coefficient 
of the host material, and R(X) is the reflectivity of the front surface. 

Gj is the generation rate of electron hole pairs through absorption of the IB, which will 
be studied later. R C v = r cv ( n P ~ n "i ) is the direct recombination between the CB and the 
VB, and rcv is the recombination coefficient. 

Rsrh is the Shockley-Read-Hall (SRH) recombination between conduction and valance 
bands, 

Rsrh = — -, — ; ^— 7 — ; r, (6) 

Tn{P +Pl) + T p (n + Tlx) 
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where r n and r p are electron and hole lifetime respectively, n\ = rijexp ( Et J^ i ) and p\ = 
n^exp { Ei uT )i n i * s the intrinsic carrier concentration of the host material, E t is the defect 
energy position, which for simplicity is set in the middle of GaAs bandgap, E^ is the intrinsic 
Fermi energy, k is the Boltzmann constant and T is the temperature of the solar cell. 

B. Electron occupation factor and net generation 

In most literatures, carrier generation and recombination via the IB are analyzed based 
on detailed balanced method, finally the results are often overestimated and deviate from 
experimental results. By adding the light generation components to the derivation process 
of SRH recombination in textbook,— the involved four processes for the generation and 
recombination can be written as:— The carrier recombination from the CB to the IB 

Rci = r C iNj(l - f)n, (7) 

the generation from the IB to the CB 

Gci^e.Njf + gcrNjf, (8) 

the recombination from the IB to the VB 

R VI = ryjNjfp, (9) 

and the generation from the VB to the IB 

G VI = e h Nj(l -f) + gviNj(l - f), (10) 

where rci and ryj are the recombination coefficients, e e and are the emitting coeffi- 
cients due to factors such as thermal excitation, g C i and gvi are light generation coefficients 
proportional to the light intensity and absorption coefficient via the IB. 

If there are no optical and electrical injection, the device is in equilibrium with a uniform 
Fermi level Efo, electron density no in the CB, hole density po in the VB and electron 
occupation / m the IB, and recombination is equal to generation, 

roiN I (l-fo)n = e e N I f , (11) 
rviN I foPo = e h N I (l-f ). (12) 
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p = rii exp ( — E/ ° ) , (14) 



Using 

n = m exp ( Ef ° kT — ) , (13) 
-E 

and 

where £7/ is the IB energy position, we can get e e and eh, 

e e = rcmi exp y~^j^J > ( 16 ) 

(£7j — -E7 \ 
— fcT — /' ^ ' 

If there is steady light injection, electrons have a stable transition from the VB to the 

CB via the IB with net transition rate Gj, 

Gi = Gci — Rci — Gyj — Ryi- (18) 
Solve the Eq. ([18]), leading tc 20 



/ = e h + 9vi + r CI n 

e e + e h + g C i + g V i + r C in + r VI p 

G _ N i [e e 9vi + e h gci + 9ci9vi - r C ir VI (pn - n 2 )] 
1 e e + e h + gci + gvi + r C in + r VI p 

This time if we set gci = gvi = 0, / and Gj can be expressed as 

/' = - e \ +rC "\ , (21) 

e e + e h + r C in + r VI p 

G > = N irci r VI (pn - nf) ^ 
1 e e + e h + r C in + r VI p' 
Gj has the form of SRH recombination, if there is no electrical injection, n = uq and p = po, 

so Eq. EH becomes a equivalent expression for Eq. [T5], 

/o = -, 1 N = • (23) 

exp ( ElzElQ j _|_ i e e + e h + r c/ n + r VI p 

The selective ohmic contact boundary conditions for Eqs. 01]) and (jSJ) are used when 
solving the drift-diffusion model, 23 J n (L) = eS n z,[n(L)—no(L)], J P (L) = — eS P L\p(L)— p (L)], 
J n (R) = —eS n R[n(R) —no(R)] and J P (R) = eS P B,[p(R) —po(R)], L (R) stands for left (right) 
boundary, Ss with subscripts are surface recombination coefficients. 
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III. SIMULATION AND RESULTS 



A. Parameters 



The simulation structure is depicted in Fig. [H the thickness of the p-region, i-region and 
n-region are 200, 500 and 500 nm respectively. The 50 nm QD layer is placed in the center of 
the i-region. The bandgap of the host material GaAs E g is 1.424 eV, the IB energy level Ej 
is formed by electron ground energy level of QD, for InAs/GaAs QD Ej = 1.124 eV.— The 
solar cell is simulated under AM 1.522 solar spectrum with the reflectivity of the front surface 
-R(A) setting to 0.1, and GaAs absorption coefficient from Ref. 2^, the carrier mobilities of 
GaAs are related to the doping, 

7200 „ , 

^ = mo+ 1 + [{ N D+ N A) / 6x wr* cm/v% (24) 

400 2 . 

^ = 32 + l + [(iV D + JV A )/1.88xlOT- 5 Cm 7 (25) 
For tqi and r VI , according to Eq. [221 the lifetimes of electron and hole are T' n = l/(Njrci) 

and Tp = l/(Nirvi), respectively. Due to Auger cooling effect,— electrons are no longer 
affected by the phonon bottleneck effect, and can easily be captured from the CB to the 
IB, so r' n is small, about 1 x 10 -12 s, Auger cooling effect is beneficial to QD lasers for high 
modulation rate, but detrimental to QD solar cells. r' p is mainly determined by electron 
radiative recombination lifetime from the IB to the VB, rj^lx 10 -9 s. So rci and ryi 
are set to 1.25 x 10 -5 cm 3 s" 1 and 1.25 x 10 -8 cm 3 s -1 , respectively. The light generation 
coefficients gci and gvi are from Ref. |20j, which are from fitting the experimental results, 
9ci — 9vi — 2.31 x 10 3 s _1 , to make a comparison, we also studied the case with g CI = 
gvi = 1 x 10 5 s -1 , which are the only parameters we assumed. Other parameters are given 
in Tab. [H 



B. N-type doping without QDs 

In order to exclude the effect of doping on the control GaAs solar cell, we first investigate 
the doping in the 50 nm layer without QDs. In Fig. [21 as the doping density increases from 
to 2 x 10 17 cm -3 , the open circuit voltage V oc increases from 0.99 to 1.01 V, while the 
short circuit current density J sc decreases from 23.67 to 23.63 mA/cm 2 , the fill factor from 
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0.82 to 0.84, and the conversion efficiency from 19.3% to 20.1%. When the doping density 
is 0, it is corresponding to the control GaAs solar cell and its data are noted as dashed line. 
The change of J sc and V oc are due to the change of total recombination Rsrh + Rev- 111 
Fig. E(a), for short circuit the total recombination between the doping layer and n-region 
increases as the doping density increases, while the total recombination decreases when the 
output voltage is IV as shown in Fig. [3](b) . In fact, the doping layer has a coulomb screen 
effect on the region between the doping layer and the n-region, the build-in electric potential 
gradually falls to the region between p-region and the doping layer, and the carriers have a 
redistribution, so the total recombination changes accordingly. 

C. Low light generation coefficients 

In the calculations below, the QDs are added to the 50 nm layer. In Fig. HI gci = 
g VI = 2.31 x 10 3 s _1 , as the doping density increases from to 2 x 10 17 cm~ 3 , V oc increases 
significantly from 0.84 to 1.0 V, for the doping density larger than 1.5 x 10 17 cm -3 , V oc is a 
little larger than that of the control GaAs solar cell 0.99 V. J sc decreases monotonically from 
23.82 to 23.65 mA/cm 2 , for the control GaAs solar cell it is 23.67 mA/cm 2 . The fill factor 
is larger than that of the control GaAs solar cell for the whole doping range, and there is a 
kink at doping density about 9 x 10 16 cm -3 , which is a little larger than Nj, it indicates the 
kink happens when the IB is fully occupied by electrons. What we most care the conversion 
efficiency increases from 17.2% to 20%, that of the control GaAs solar cell is 19.3%. 

D. High light generation coefficients 

To better understand the effects of the doping on current-voltage characteristic, we cal- 
culated another set of data for g C i and gvi, gci — 9vi = 1 x 10 5 s _1 , corresponding to high 
photon absorption or wide range absorption spectrum. The results are shown in Fig. [51 V oc 
has a large increase, which is almost the same to that for low generation coefficients, while 
J sc has a large increase with maximum value 30.07 mA/cm 2 , it decrease to 23.66 mA/cm 2 
as the doping density increases. The fill factor has a large fluctuation, its value is as low as 
0.74 at doping density about 9 x 10 16 cm~ 3 , at which doping density there is also a kink. The 
conversion efficiency has a large increase compared to the result in Fig. 11(d), its maximum 
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value is 21.6% when there is no doping, it decreases to the smallest value 18.8% at about 
doping density 8 x 10 16 cm -3 , which is equal to Nj, then increases to 20% as doping density 
continues to increase. 

Comparing the results above, we can draw the conclusion that when the doping density is 
large, such as doping density greater than 1 x 10 17 cm -3 , the increase of conversion efficiency 
is mainly attributed to the doping effect on the control GaAs solar cell without QDs. This 
is due to the IB is fully filled with electrons for the high doping density, and there are no 
empty states to accept electrons for tansition from the VB to the IB, and the role of QDs 
are suppressed. 

For the QD IBSC with high generation coefficients, Fig. [6](a) and Fig. [6](b) show n and p 
distribution in the host material GaAs for short circuit and output voltage 1 V, respectively, 
and the electron occupation factor / in the IB is depicted in Fig. [7J When there is no 
doping, / is small on short circuit condition, —f'Nj is small in Eq. [21 so the small amount of 
electrons in the IB have a little coulomb effect on n and p. When the output voltage is 1 V 
without doping, / is about 0.41, the electrons in the IB has a large coulomb effect, and lead 
to a concave in n and a convex in p as seen in Fig. Mb). When doping density increases, 
under short circuit n in the region between the QD layer and n-region increases evidently, 
for the output voltage IV, n increases and p decreases at the position of the QD layer, and 
there is a convex in n and concave in p for large doping density. 

E. Analysis and discussion 

It is easy to understand the change of J sc and V oc as doping density increases by simply 
analyzing Gi. As shown in Fig. El n and p change as the QD layer's position, the doping 
density N dope and output voltage V op . According to Eq. [201 Gi is a function of n and p, 
considering the QD layer is fixed in the middle of the i-region in our study, so Gj can be 
written as G/(A dope , V op ). 

On short circuit condition, if there is no doping, the electron and hole carrier concentra- 
tion are relatively small in the middle of the i-region as shown in Fig. E(a), so we can neglect 
the items about n and p in Eq. ( 1201) . therefore, 

g, W». = o, v, = o) » Nl(e ' 9v ' + + gcm,) . (26) 

e e + e h + gci + gvi 
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Because e e 3> e^, e e ^> gci and e e 3> gvi, Gi(Nd ope = 0, V op = 0) can be further simplified, 
Gi(Ndo pe = 0, V op = 0) ~ Nigvi, this means Gj is fully determined by gvi on short circuit 
condition without doping, compared to the control GaAs solar cell its increase AJ SC is 
0.15 niA cm -2 in Fig. HJb) and 6.4 mA cm -2 in Fig. 0(b), fully conformed to the relation 

AJ SC « eGjWQo = eNjgyjWQD, (27) 

where u>q,d is the thickness of the QD layer 50 nm. When considering the doping, n at the 
position of the QD layer increases, so n can not be neglect, 

GHiW.^O)* *''^'-^™' . (28) 

e e + r C in 

Apparently, Gi decreases monotonically as the doping density increases, so J sc decreases 
monotonically. 

On open circuit condition, according to Fig. Mb), n and p are large, considering e e 3> e^, 
e e > gci and e e > g VI , so 

G,(^.V^ = y.)»- f r " r f' . (29) 

e e + rc/n + r y/ p 

Due to the increase of n and tqi 3> ry/, gets smaller, that is the recombination via IB 
is reduced, so V oc increases. 

By now literatures all point out that the IB should be half-filled.— Jii22i2£ But in our 
simulation model, there is no such relation between the conversion efficiency in Fig. [5^d) 
and its corresponding electron occupation factor in Fig. [71 With the increase of the doping, 
the predicted increase of light current does't appear in our model. Because the maximum 
short circuit current density is determined by Eq. (27J and can't be increased more by doping. 
Additionally, by comparing Eqs. [23] and [T9J although the IB can be half-filled by doping, 
fo = 0.5, / is no long 0.5 when the solar cell works at the maximum power poit. 

It is noteworthy that the conversion efficiency are also affected by the fill factor, in Fig. [5] 
when the doping density changes from to 5 x 10 16 cm -3 , the short circuit current density 
almost has no change, but the conversion efficiency decreases due to the decrease of the fill 
factor. 



IV. CONCLUSIONS 



In conclusion, the current-voltage characteristic of the IBSC is affected by the n-type 
doping density in the QD layer. The open circuit voltage, short current density, fill factor 
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and conversion efficiency all vary with the doping density. As the n-type doping density 
increases, the open circuit voltage increases, while the short circuit current density decreases, 
and the conversion efficiency tends to close to that of the control solar cell. In one word, the 
n-type doping tends to suppress the role of QDs, whether the QDs originally play positive 
or negative role. This prediction has been observed in our recent experiment, the details 
about growth, fabrication and experimental data will be published later. 
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TABLE I. Simulation parameters 
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P i n 


200 500 500 nm 

FIG. 1. (Color online) The simulation structure, the 50 nm QD layer is placed in the middle of 
the i-region. 
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FIG. 2. (Color online) Doping without QDs, the current-voltage characteristic changes as the 
doping density increases, (a) Open circuit voltage, (b) Short circuit current density, (c) Fill 
factor, (d) Conversion efficiency. Data for the control GaAs solar cell without QD layer are 
plotted as dashed lines. 
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FIG. 3. (Color online) Doping without QDs, the total recombination Rsrh + Rev changes as 
different doping density, (a) On short circuit condition, (b) When the output voltage is 1 V. Data 
for the control GaAs solar cell without QD layer and doping are plotted as dashed line. 
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FIG. 4. (Color online) For gci = gvi = 2.31 x 10 3 s _1 , the current-voltage characteristic changes 
as the doping density increases, (a) Open circuit voltage, (b) Short circuit current density, (c) 
Fill factor, (d) Conversion efficiency. Data for the control GaAs solar cell without QD layer are 
plotted as dashed lines. 
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FIG. 5. (Color online) For gci = gvi = 1 X 10 5 s , the current- voltage characteristic changes 
as the doping density increases, (a) Open circuit voltage, (b) Short circuit current density, (c) 
Fill factor, (d) Conversion efficiency Data for the control GaAs solar cell without QD layer are 
plotted as dashed lines. 
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FIG. 6. (Color online) For gcj = gvi = 1 X 10 5 s _1 , the carrier concentration of electron and hole 
in the GaAs host material for the doping density 0, 5 x 10 16 , 1 x 10 17 , 1.5 x 10 17 and 2 x 10 17 cm~ 3 . 
(a) On short circuit condition (b) When output voltage is 1 V. The dashed arrow indicates the 
change direction as doping density increases. 
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FIG. 7. (Color online) For gci = gvi = 1 x 10 5 s -1 , the electron occupation factor / in the IB vs. 
output voltage V op with different doping density. The dashed arrow indicates the change direction 
as doping density increases. 
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